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COMPRESSIBLE FLOW ABOUT SYMMETRICAL JOUKOWSKI PROFILES 

By Gael Kaplan 



SUMMARY 

The method of Poggi is employed for the determination 
of the effects of compressibility upon the flow past an 
obstacle. A general expression for the velocity increment 
due to compressibility is obtained. This general result 
holds whatever the shape of the obstacle; but, in order to 
obtain the complete solution, it is necessary to know a cer- 
tain Fourier expansion of the square of the velocity of 
flow past the obstacle. An application is made to the case 
of a symmetrical Joukowski profile with a sharp trailing 
edge, fixed in a stream of velocity Vq at an arbitrary angle 
of attack and with the circulation determined by the Kutta 
condition. The results are obtained in a closed form and 
are exact insofar as the second approximation to the com- 
pressible flow is concerned, the first approximation being 
the result for the corresponding incompressible flow. For- 
mulas for the lift and moment analogous to the Blasius 
formulas in incompressible flow are developed and are 
applied to thin symmetrical Joukowski profiles for small 
angles of attack. 

Since actual experimental data for Joukowski profiles 
are lacking, the theoretical results are applied to a thin 
and a thick profile at zero angle of attack, and the velocity 
and pressure distributions are calculated and compared 
with those for the corresponding incompressible cases. The 
critical values for the ratio of the stream velocity v 0 to the 
velocity of sound in the stream c Q , corresponding to the 
attainment of the local velocity of sound c by the fluid on 
the surface of the airfoils, are also obtained. 

INTRODUCTION 

When a compressible fluid streams past a fixed body 
with a velocity small enough so that nowhere in the 
fluid is the local velocity of sound exceeded, the flow 
may be represented by a velocity potential. The effect 
of compressibility is to distort the streamline picture 
associated with the corresponding incompressible flow. 
This distortion has been calculated by Janzen (refer- 
ence 1) and Rayleigh (reference 2) for circular cylinders 
and spheres and recently for elliptical cylinders by 
Hooker (reference 3). The methods used by these 
authors, however, are not feasible for the determination 
of the flow about obstacles other than the simple ones 
mentioned. On the other hand, a method introduced 



by Poggi (reference 4) may be used in determining the 
flow about shapes resembling airfoil profiles. 

The method due to Poggi is as follows: When the 
fluid is compressible, the equation of continuity may be 
written as 

bv x bv, = 1 Dp m 
bx^by pDt w 
where the symbol DfDt denotes, as usual, the operator 

v *> v >/> nu *d velocity components; 

and p, the variable density of the fluid. 

This divergence will introduce extra terms in the 
expressions for the velocity components, the divergence 
at an element dxdy being equivalent to a simple source 
_lDp 
P Dt 



of strength 



2s- 



-dxdy. Poggi thus replaces the com- 



pressible flow by an incompressible flow due to a suitable 
distribution of sources throughout the region of flow. 

If the motion of the fluid is steady, then the equation 
of continuity and Euler's differential equations of 
motion become: 



and 



bv x [ bo r Pjdp Vjbp 

bx + by~ pdz pby 

„ dflr . bv x _ Ibp 
v 'bx +v 'by~ pbx 

bv, , bv r lbp 



(2) 



bx^ "by pby 

Assuming the pressure p to be a function of the density 
p only and introducing the local velocity of sound 

t 

,bv x . fbv* , bo,\ , ,bv r 



equations (2) yield the following: 

bv x bv r _ Ij~ 
bx + by~cH 



or if 
then 



•o t =v x -\-v* 



bv x , bv. 



bx by 

If it is further assumed that the motion of the fluid is 
irrotational, then 



bv £ _bv x=Q 
bx by 
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and a velocity potential 4> may be introduced, where 
„ _ „ _ a * 

The strength of the source at a point (x, y), given by 
the expression on the right-hand side of equation (3), 
then becomes 

4vc\bx dx byiyyj 

Suppose now that (£, rj) and (x, y) are the rectangular 
coordinates of points in the f and z planes, respectively, 
and furthermore that these two planes are conformally 
related, that is 

where ?=£+ii?, z=x+iy. Let the £ plane be the 
plane of the profile and the z plane, the plane of the 
circle into which the profile is mapped by the foregoing 
conformal transformation. It is well known that, at a 
pair of corresponding points at which £ and z possess no 
singularities, a source at one such point corresponds to 




Figure 1.— Image of a simple source with regard to a circle. 

a source of equal strength at the other. It follows then 
that at corresponding elements 



1 



4xc : 



4:Tc\dx dx dy dy 



(4) 



where, in the expression on the right-band side, <t> is the 
velocity potential in the z plane while v is the magnitude 
of the velocity in the £ plane. 

In polar coordinates (r, 0) the strength of a'source at 
an element dxdy of the z plane is 



> 



or, introducing a new variable X=— (where if is the 
radius of the circle into which the profile is mapped), and 

»r=— |j~> p»= — ~ this expression becomes 

1 / do 1 v, W\K JSJa 
■£A V ' bX~X W)\ d}dd W 

With the source distribution known in the plane of 
the circle and given by equation (5), the induced tan- 
gential, velocity at the circular boundary may be cal- 
culated. 

Thus consider a unit source located at a point Q of 
the z plane. In the presence of a circular boundary of 
radius B, the velocity induced at any point P external 
tor or on the boundary is given by 



dw = _/_l _1 1\ 

~d~z \z—Zq z—z a z) 



where S is the point inverse to Q in the circle. (See 
fig. 1.) Since the normal velocity at the boundary is 
zero, the velocity there is wholly tangential and is 
given by 



C 



dw\ 2X sin (9— S) 

f[z~/ un fl[l-2Xcos (fl-8) + X"] 

where z=i?c ,{ ; ZQ=re u ; and z a ——t* t 



(G) 



Hence, the total velocity induced at any point of the 
circular boundary by the system of sources given by 
equation (5) is 



V 2«%Vo Jo 1-2X 



dv 1 Vf So* 

-2 X & cos(g-^ + X^ 8itl («-«^C7) 

The justification for replacing c by c 0 in equation (7) may be 
shown in the following way. From the Euler equations of 
motion (2) and the condition for irrotational motion, It follows 
that 

Then when adiabatio conditions prevail so that the relation 
between p and p Is 

it follows by integration that 



where the zero subscripts denote the corresponding magnitudes 
in the undisturbed stream. From the foregoing equation it is 
seen that e has a maximum value at the stagnation point where 
»=0, that is 
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Furthermore, as the streamline corresponding to the boundary 
of the obstacle is traversed, a point is reached where, for a definite 
value of the ratio pj/co, the velocity of the fluid equals that of 
the local velocity of sound. This critical velocity is obtained 
from equation (8) by putting c=r and solving for v. Thus 



(9) 



For example, let ro/co=0.75. Then with 7=1-408 (for air) 

C,,a*=1.056 Cfc 



and 



C»f»=»cr£t= 0.062 Co 



Away from the obstacle, p approaches Vq and c approaches <h. 
Thus it is seen that the variation of the local velocity of sound 
from c, is, in general, small enough to permit replacing c by Co, 
at least to a first approximation. 

Equation (7) is a functional equation for the fluid 
velocity v and may be solved by a method of succes- 
sive substitutions. The procedure, due to Poggi, is to 
substitute for v T) v>, and & values pertaining to the 
corresponding incompressible flow and thus obtain a 
first approximation to the sink-source distribution in 
the plane of the circle. The method thus considers 
the incompressible flow to be the first approximation 
to the compressible flow. The second approximation 
is then obtained by superposing on the incompressible 
flow the effect of the sink -source distribution as given 
by equation (7); that is 

i , «»p=y<«« m j+A» (10) 

GENERAL DEVELOPMENTS 

Before equation (7) is applied to any particular case, 
it is expedient to consider it first in a general way. 
Thus, suppose that c* can be developed in a Fourier 
series so that 

r s 1 " 

-j=o«o+S(ff» cos 710+6, sin Tiff) (11) 

where the a nt b t are functions of X and also contain the 
parameters of the shape. 

Then 

dX 



^=Po ! [^o'+S Co.' cos nfl+6/ sin 7i0)] 



be 



=ni a y,7if6. cos nd — ce* sin nd] 



Then the potential is 



and 



Po 



dz 



z 1 



(12) 



Then 



where the primes denote differentiation with regard 
to X. Expressions for v r and r» are obtained from the 
complex potential of the flow past a circular cylinder 
of radius S, with the circulation determined by the 
Kutta condition. Let the stream of velocity p 0 make 
an angle a with the negative direction of the x axis. 



tV=— X 2 ) cos (fl+c£)=!? 0 (ci cos 0+Cj sin 0) 

and 

ti=v a (l+\ T ) sin (0+a)-f-2cbXsin a 
=vj^+dt cos 8+d 2 sin 0^ 

where 

<?!= — (1— X s ) cos a 
ci==(l— X s ) sin a 

"2"=2X sm a. 

d 1 ={l+\ r ) sin a 
^2=(1 + X 2 ) cos a 

Therefore 

p,^=-| EbWfo cos 8+Ci sin 0) 

+5*o 8 S{a«'c 1 [cos (n-l)0+cos (n+l)0] 
* n-l 

+6.'c 1 [sin (re+l)0+sin (n— 1)0] 
+a/cjsin. (n+ 1)0— sin (n— 1)01 
+b x 'c 2 [co3 (n—l)0— cos (n+l)0]} 



and 



v, bv 2 



— Wilcos (n— l)0+cos (n,+l)0] 

+ 2«^2[C0S (?l — 1)0— COS (71+1)0] 

— &„<f 2 [sin (n+l)0— sin (ti— 1)0] 
+a*d a sin nd—b^o cos nd} 

The following definite integrals will be found useful in 
evaluating equation (7) : 

ran (0-5) fl^ fl _° if ?i=0 

1-2 X cos (0-5) + X iSm 710 aa ~ a-X"- 1 cos nt> if 1 



X 



* sin (0—5) fl^ fl _° if n.=0 

o l-2Xcos (e-5) + x* COS7lffcZ& - -xX-'smnSifngl 

dr 2 

Then substitutmg the foregoing expressions for zv - ^ 
and y -g^ into equation (7) and integrating with re- 
gard to 0, it follows, after replacing the derivatives 
a/, W by a„, b x by means of partial integrations, that 
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r 

CO fl 

+ (&s)x_o 009 ( s_ a ) — S2n cos (riS+a) X*6 s _idX 
n-3 «/0 

+S2»cos (nl—a) Px n - 2 £> n+1 a*X 

n-l Jo 



+£>&sin (7i5+a) fVc-idX 

n-l Jo 

-S2?3,sirL(7i5-a) Px*- 2 a„ + i<^ 
n=I Jo 

<» /*l 

+sin oiT^An sin nSj^X"~ I J ll dX 



+sin a^Zin cos nfi f X"~ 1 a B d\ \ 

n-l Jo J 



(13) 



f plane, and the circle of radius R with center at the 
origin 0 of the z plane into a Joukowski profile with a 
sharp trailing edge in the f plane. (See fig. 2.) The 
distance 00' is denoted by ae and e is therefore a measure 
of the thickness of the profile. Since the profile has a 
sharp trailing edge, the two circles touch at the cor- 
responding point (—a, 0). The relation between the 
z, z' planes is 

z'=ea+z 



z, z' planes 




If w denotes the complex potential of the incompres- 
sible flow in the £ plane, then the complex velocity is 
given by 

dw = dw dz dz' 
d{ lzdT'~dJ 

where .. . 



$ plone 



x,x' 




Figure 2.— Transformation of a symmetrloal Joukowski profile Into a circular contour, 

dz 



It is to be especially noted that this expression for At? 
is perfectly general and independent of the type of pro- 
file considered. All that is required for a complete 
solution of equation (13) is a knowledge of the Fourier 
development for v 2 lv 0 2 . 

APPLICATION TO SYMMETRICAL JOUKOWSKI 
PROFILES 

Suppose now that the boundary in the f plane is a 
symmetrical Joukowski profile with a sharp trailing 
edge. The Joukowski transformation given by 



•. I r a ~ 



(14) 



maps the circle of radius a with center at the origin 0' of 
the z' plane into the line segment (—2a, 0; 2a, 0) in the 



gp-l and -fi= iz > +a)iz ,_ a y 



According to equation (12) 



Hence 



But 



Therefore 



dw_ (z+R)(ze (a -Re-'°) 

dw_ (z+R)(ze t «-Re- t °)z ,t 
~37~ Vo z t {z , +a)(z , -a) 

z'+a=z+a(l+e)=z+R 
z / — a=z— <z(l — e) 



dw (z+ea) s (ze { "— Re~ la ) 
d^ V ° z*[z-aQ.-6)] 
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Introducing 



. R ea , , a(l — e) , 
X= P E =h ' and— 5—=* 



it follows that 

i (l+2AXcosfl+A 1 \ i ) 2 [l-2X cos (fl+2a)+X*I 
m 1— 2JtX cos 0+FX 2 

It is now required to obtain the Fourier series for 
t^/Po*. Thus, making use of the following developments: 

sin B A,,.., , . 

1-2*X eort+Ptf -a^ ^ 719 

l-2ifcX cos «+t"X»=T=3^ i+2 §* X) * C0S nd ] 



it follows that 

£-<l+2AX cos g+A'X ^ + ^T^x' 
+2X sin 2aS(fr*)" _1 sin 



cos 2a 



n-l 



G 



fc(l + X 2 )-(l+PX 2 ) cos 2a - 

+ 2 *(i-*»x") § (AX) cosn * 

Reducing this expression to the form of a Fourier 
aeries, it turns out that, for n>2 



fl,=2*-« (A +*)* (1 +A*X*)*X" 

r „ , £(l+X 2 )-2 cos2al 

L cos2a:+ i=h? — J 

6.=2i"- s (/i+A:) t (l+M-X t ) i X» sin 2a 



(15) 



For later use it will be convenient to introduce the 
following notation for n^O: 



a„ l =2*"-»(A+*) 1 (l+A*X 1 )V cos 2a 
a, s =2*»-»(A+i) 4 (l +MX*)»X ^' (1 + ^~ 2 X , C0S 2 " 



6, =2&*- 3 (A+£) 2 (l+AjtX 2 ) 2 X' sin 2a 
Also 

^a 0 =AD+ (5+ Ck\)Fk\ 

a,=2BD+ CA+BiX+Cfc^iftX 
a»=2CD+ CB+iliX+B^+GPX^fJfeX 
fc^-Sfcl+SA-X+C^X 2 ) 

where 
4=1+4A S X 2 +A«X 4 
S=2AX(l+A 2 X i ) 
C=A 2 X 2 



(16) 



(17) 



rt__ l + (l— 2k cos 2a) X 2 
U l-PX 3 

.E=2Xsin2a 

r ,_2(A--cos 2a)+2fr(l-fc cos 2a) X 2 
£(1-PX 2 ) 

It is seen that 

(ao)x-o=2 and (ai)\-o=(6a)x_o=0 
Equation (13) then becomes 
A»=^£-2 sin(5+a)+2 cos (?-a) J^bt-bJdX 

+4 cos (25— a) j\(bt—l> 3 )d\ 

—2 cos (5+a) Pfo-SOrfX 
Jo 

—4 cos (25+ a) Px(5»— &,)«*X 
Jo 

-4 cos (25+a)J^ 1 X 2 (& 1 -5 I )rfX 
—6 cos (35+ a) j\ s (b 2 -l s )d\ 
+2 cos (5+a)J^XMX 
+2 sin (5+a)J^(<n— «i)<£X 
+4 sin (25+a) J\(a g — a 4 )<ZX 
—2 sin (5— a) (a x — a{)d\ 

—4 sin (25— a) \\(a t — a^rfX 
Jo • 

+2 sin {8+a)j\(a 0 —a a )d\ 
+4 sin (25+a) J\ J (a 1 -a I )dX 
+6 sin (SS+aJ^X^Cj-aO^X 
+2 cos (5-a)J^(6,-&,)rfX 
—2 sin (5— 0)^^(03— a a )cJX 
-22ft cos (n5+a) J^X-J^+X- 1 !.)^ 



(18) 



+S2n cos 



+S2n sin 



n-l 



l 6^ 1 +X- I 6»)rfX 
^ 1 +X»~ 1 a,)dX 



-S2ft sin (n5-a) f V^a.+i+X-'a^xl 
*~i Jo J 
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S2nsin (ni-a) f\x»- , a II+ , , + X»- 1 a 11 , )dX 

n-l Jn 



202 
where 

Oo-^-p {(&+*)•(«» 2«-*)+Jfe»+[ft(A«+jfe") 

-*(l+A it )(/i+A:) s +2M(A+*) 2 cos 2a 
+2A 2 F cos 2a]\*+[-h 1 k(h+ky 
+h l k+h i k 2 (h+k) 2 cos 2a]X 4 } 

2A*X l + X 2 (l-2/fc cos 2a) 



1-FX 2 



- 2A S 
Os— ai=-r- X s cos 2a 



(19) 



&1-&1 2X ™ 2a ft(H-2*) +2A ! /:X 2 ] 

63—63= — ~£-X 2 sin 2a 

6 0 =p (A+&) 2 (1+A*X 2 ) 2 gin 2« 



It is a great simplification to replace a„ by a^+a^ in 
the last two integrals of equation (18) before perform- 
ing the integrations. Then 

-it^n cos (nS+a) f\x»6 II _ I -|-X"- I 6 Il )dX 

s-l JO 

+ S2ntrin (nS+a) f '(XX+i'+X'-V^X 
n-l Jo 



= -(l+A-)(/ t +A:) 1 24n^ tt - < sin (n5-a) 

B-l 

ftt+Ajfcx 2 ) 2 * 1 - 1 ^ 

+ (1-2* cos 2a+P) (1+jfc) (h+k)* 
2j4n£ sin (nS— a) i- — rfx 

n-l Jo 1 — ft X 

Consider. now the integrals 

Ji=Sn*"ain (n«+|8) f 1 (l+/iiX 1 ) 2 X 2 » , -WX 

11-1 Jo 



(21b) 



and 



Then 



n-l Jo 



and 



= (l+fc)(A+&) 2 2>nfc*- 4 sin (na-o) 

n-l 

f I (l+Ai-X 2 ) 2 X 2 »- 1 ciX 
Jo 

S2»cos InS-a) f(X fl - , 6,, +1 -r-X n - 1 6,,)<zX 
n=i Jo, 

-f$n sin(n5-«) f\\ m -Hh+i l +\' F W)d\ 

n-l Jo 

=-(l+fe)(A-f-^:) 2 S4nJt' , - ,l sin (n«-3«) 

n-l 

J o I (i+Mx 2 ) 2 x 2 »- I dx 



(20) 



Also 



n-l Jo 

= - (1 +k) (A+ifc) 2 S4nt"- 5 sin {nS+a) 

n-l 

J Q 1 (H-AitX 2 ) 2 X 2 ' , - 1 (fX 

+ (1 -2k cos 2a+£ 2 ) (1 +k) (A+£) 2 

S4nfc- 5 sin (««+«) P^i^^^x 

n-l Jo 1 — * A 



(21a) 



«)i+tJi=fl*S +/ii-X 1 ) , nl-''X 2 ' 1 - 1 c < "*dX 

n-ljo 

Let s=l-X 2 . Then 
r 1 .r-^r ^^+^V - ^'^T ^f- 1 (1+**-*)' 

~e"(l+A<r«) s +2A(H-A<r«) log 

Hence 

r Kl+Afr). 2 an (g+fl-fr sin ft 
Jl ~ 2 1-2* cos 8+ A 1 



— A^cos (5— 



„, . _, a: sm 0 

/S) tan ^ — j , 

1— £cos8 



tsinS 

-t-^sin (5-/5) log (1-2A cos 5-j-A- s )-£ sin 0 



-}?\k sin (5-/3) + cos (28-/8) tan-' ^j^ g 

+|sin (25- 0) log (1-2* cos 3+k 1 ) —k 1 sin jsj 

Replacing /3 by a, -a, or —3a, the corresponding integrals 
in equations (20) and (21) are obtained. 
Consider finally the integrals 



: - f I (l+AJtX ! ) 2 X s «- 1 ,, 
ii=Z/^ cos nS I izrj^i ^ x 



and 



Then 



r ^ 7» • , Pd+AitX^X 2 *- 1 
am n*J o ^jr^ 



^X 
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Again let s=k\ 2 , 

, . ., e a C k (1+A s)» , 

fi+u*— 2j 0 (l-4 8 )(l-se«)^ = 

-!3? logCi - i2) 



(A+4)(A+f) 
+ 2(4-e«) 2 g ( 
_4(A+e«) i _ 
2(l-4e*)(4-e <i ) 

Therefore 

il_ 24 (1-24 cos5+4 i ) jiogei_M 
A (A4— 1) cos 5 +4— A cos 25 



■ke a ) 



1-24 cos 5+4* 
A (A4— 1) sin 5— A sin 25 



tan' 



log (1-24^3 5+4*) 

_ t 4 sin 5 
1—4 cos 3 



1—2* cos 5+4* 

A+4 [A4*-24+ (1 -24A+4*) cos 5+A cos 25] 
4 (1-24 cos 5+P)* 

log (1-24 cos 5+4 2 ) 

A+4 (1 — 24A— 4*) sin 5+A sin 25 x A sin 5 



(1-24 cosS+4 2 )* 



1— A- cos 5 



+5 



£ (l+A*) cos5+2A 



2 1-24 cos 5+4 1 



and 

r_ (A+4)' (1-4*) sin 5 , 

J * 2F~ (l-24cosS+4*)* logC1_F) 

. A (A4— 1) sin 5— A sin 25, nT . . ... 

+4 1-24 cos 5+4* l0g (1 ~ 2£ C0S 5+ ^ } 



+5 



A (A4— 1) cos 5+4— A cos 25 ^ _i 4 s in 5 



1-24 cos 5+4* 



1—4 cos 5 



A+4 hk i -2k+ (1-2A4+4*) cos 5+A cos 25 
2 (1—24 cos 5+4*)* 

k sin 5 



tan" 



1—4 cos 5 



A+4 (1— 2A4-4*) sin S+A sin 25 
4 (1—24 cos 5+4*)* 

log (1—24 cos 5+4*) 

4 (1-A*) sin 5 
" f "2 1—24 cos 5+4 1 



From equations (18), (19), and (20) and the integrals 
Ji, I u and Ji, it follows that 



(9+5A*)(A+4)' 
3F 



3A* 5 A* 2A(A+24) 



4 

(3+3A4+A*4*)(A+4)» 



4 s 



T*~ 
2A 1 
"T 



■] 



16 A* 



sin (5+«)+[§ 

+ 2A(A+24) j sin (s+3a) _-» ^ (2J+a) 

. [~10A* . 2A(A+24)"1 . ,„.,, w 2i'. , OI . _ . 

+ L3 £ — — J sm ( 2s + 3a )+T sm ^ 

. ^ ^ r* \ , 8(1+4) (A+4)' r , . 
5111 (5-t*)+ ' ^ Ji(— g) 



,4(l+4KA+4)^ c _ 3ffl) . 



4(l+4)(A+4)' rr > 
p W 



+ A(i+ky(h+ky {1 _ 2k cos 2a+p) 7i sina 
+ 4(i-FKA+4) i (1 _ 2fc cog 2a+AJ) ^ ^ 1 (22) 

where J 2 (—a), for example, means that in the expres- 
sion for Ji (/3), —a has been substituted for j8; and 

' ^ere ^ no difficulty in evaluating Ap/t>„ for 

any value of A except A=0. For A=0, 4=1, however, 
the JoukowsM profile degenerates into a line segment 
and 



^=>M sin «Tcos(5+2a)+2 sin 4*9+*)+*** 
v„ L 1— cos 5 

log2(l-cos5)- 1 i^i)Il 
1 — cos Sizinj 



2 sin 8 a 
1— cos 5 



(23) 



The last term in this expression contains the divergent 
" 1 

series S - ' which approaches infinity like T.ttti log n. 

0=1 n B-»oo 

This infinite term shows that streamline flow cannot be 
maintained about a straight-line profile except for the 
trivial case of zero angle of attack. 

CORRECTION FOR THE CIRCULATION 

It is noted that the expression for the complex 
velocity about the circular profile given by equation 
(12) was obtained with the circulation fixed by the 
Kutta condition. When to, representing the effect of 
compressibility, is added to the incompressible velocity 
obtained from equation (12) to yield the compressible 
velocity, the Kutta condition no longer holds. In order 
to restore the Kutta condition, an additional circulation 
Ar is added to the incompressible one. Thus, the 
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Telocity at the boundary of the circular profile is given 

by 

o • I \ I r 0 +Ar . ad 

-=2 em ( S+ «)+^-+- 



(24) 



where r 0 +Ar=r and Y 0 =4irRvo sin a, the circulation 
in the incompressible flow. The Kutta condition, i. e. , 



0, thus serves to evaluate Ar. For S=v the 



expressions for J 2 (/3), I u and 7 S simphfy considerably. 
Thus 

(Jt(g))»-r=- f 2 (1 (1 ^' - l A&og 

+A ! [A-|A 2 -log (l+i)])siii/S 

(/i),-.=|log (l+^+gglog 
(Js).-.=0 

The relation between the velocities in the plane of the 
circle and the plane of the profile is given by 

dw_dw dz dz' 

~d~f~dl dz''di 

a 2 

where z f =z+ae and i"=s'+p- 
Then 



dz 



dz' 
dz 1 



= 1 and 

_ l+2&Xco3,3 +A*\* 
Vl+2Xcos5+X 2 Vl- 



-2&X cos 5+^ 
It follows that on the profile where X=l 

Vp__ 1+2A cos v e 



(25) 



«o ' V2 (1 + cos 5) VI —2k cos « +£ J »o 

where t> P is the velocity on the profile corresponding to 
v e on the circle. 

When the profile is assumed to be thin so that only 
the first power of A is retained and the angle of attack is 
small so that cos ae*l and sin a^a, then the Kutta 

AT 

condition leads to the following expression for 



2^=^(1 +h)a 



It then follows that 



r , ,1+.A 



(26) 



(27) 



This value for the ratio T/r 0 corroborates Glauert's 
result (reference 5) 

.,1 , 

l l+2^+ • * ' 

when the profile is very thin, i. e., when h is negligible 
in comparison with unity. 

Since the rigorous expressions are available, it may 
be interesting to compare the approximate result 



given by equation (26) with the exact result obtained 
from equation (22). Thus, for a very thin profilo 



more 



defined by e=0.01^A=— jji A-=^^ and for the 
conventional profiles defined by e=0.05^A=^i £=^y) 
and €=0.10^A=-yi k—j^j, at angles of attack a=10° 
and a=5°, the following table presents the results: 





•-0.01 


1-0.05 


t-a ic 




Ar 

inRtt 

(approxi- 
mate) 


Ar 

2t&>o 
(exact) 


Ar 

■hrRct 

(approxi- 
mate) 


Ar 
(exact) 


Ar 

(approxi- 
mate) 


Ar 

2rftl 

(exact) 


10 

s 


0. 17«M> 
.OSglOja 


.0921S/I 


0.15281m 
.09U% 


0.20281m 
.091%! 


aiwai/t 
.oosM/i 


.03557,1 



It is to be noticed that for a= 10° the exact evaluation 
of Ar/2xi?ito yields a greater value for e=0.01 than for 
e=0.05, a fact not given by the approximate equation 
(26). This reversal appears, in general, for larger values 
of e as the angle of attack increases; e. g., for «=0.05 
at a=20°. This feature of the exact expression for the 
additional lift has no practical significance, however, 
insofar as the lift is concerned, since the appropriate 
combination of e and a showing this reversal is outside 
the practical range. 

In the calculation of the local velocities and pressures 
on the surface of the airfoil, the rigorous expressions for 
Ap/co are to be used. The rigorous derivation, however, 
of the total integrated lift and moment ou the airfoil 
involves great mathematical difficulties. A simplified 
form fox At)/» 0 may, however, bo obtained for a thin 
Joukowski profile at small angles of attack. Its use in 
integrating for the lift yields, as will be shown later, 
the expected result that 

Lift =p 0 for 

where V=T<(l + 1+h 

This result justifies the use of the simplified form of 
Av/v 0 in calculating the lift, but its use in integrating 
for the moment, although reasonable, is somewhat 
uncertain. 

If, then, only the first power of A is retained and the 
angle of attack taken small enough so that cos a=*l and 
sin aas«, it follows from equation (22) that 

-^{[cos 5+A(4+8 cos 5+3 cos 25)]a+A(sin 5+sin 25) } 

The expression for cjvo, replacing Ar/2T2?iio by tho 
value given by equation (26), then becomes: 

- e =2sin5+2«(l+cos5) 

+/t[(l+cos 5)a+A(5+8 cos S+3 cos 2S)a 
+A(sin6+sin25)] (28) 
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CALCULATION OF THE PRESSURE AND LIFT ON THE 
AIRFOIL 

According to equation (8) 
Then from the adiabatic equation of state 
and the definition of the local velocity of sound c 



it follows that 
Therefore 



[dp 



c? (po) Co) 7 



(29) 



Expanding the right-hand side of the foregoing equa- 

tion according to powers of — (=/«) it follows that 

Co 

P=Po+^o 2 (l-^)-r~A,ro i ^(l-J)+ • • • (30) 

The pressure distribution may be calculated by means 
of equation (29) together with the values for c/r 0 obtained 
from equations (22), (24), and (25). Equation (30) 
will be used in obtaining the total lift and moment on 
the airfoil. 

Since the profile is a streamline, the normal velocity 
— d«£/dn=0 and, accordingly, if s denotes the length 
along it, then Bernoulli's equation may be written 

J ,-con 3 tan t -i R ( 1+ I,)(0 + i^)' + . . . 

Let n denote the inward-drawn normal to the con- 
tour. Then from figure 2, it is seen that p cos (£,7i) 
and p cos (ij,n) are the components of the pressure along 
the £ and 17 axes, respectively. Accordingly, the force 
on the airfoil is given by 

F=Pe-iP,=-f (l+|,)^J[cos 

-i cos ( V ,n)]d8+1^ "£(57) W &n) 

—i cos (ri,ri)\ds 

where the profile G is traversed in the counterclockwise 
positive sense. On the other hand 

d%=ds cos (i?,n), dii=ds cos [t— (£,7i)]=— ds cos (f,n) 
and therefore 



Now, by definition, 
dw 



and, since the velocity normal to the profile equals zero, 
it follows that 

g_^ (S ,»)|i_^<,,„) 



Therefore 

But 
or 

Hence 
and 



dw 5<fr d£—idri 

df~bs ds 

ds i =d?+dr?= (dt+idr,) {di-idrj) 

d£—idri=^ 

i>4> = dw df 
bs df 5s 



Now 



0j0 =^ where f is the conjugate of f 
Therefore 

Referring to the plane of the circle of radius J? 

-^L(g)'(l)'l^ _ « 

where z=e" and 5 is the polar angle of the circle of 
radius R. 

Since by definition 

it follows from equation (28) that 

^=mw-«=Bv c (a 1 z+a 0 +^+^+'%) (33) 

where 

a„ =(l+^)+i«(l+f+^) 
o_i=ta(2 + ^+5A/i) 

a -* — ( i +^)+K l+ f +4Ajct ) 

hfl . 3 .r 

£t-8=— 2+2 P a 
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But 



Also from the Joukowski transformation 



and the relation 

it follows that 

dz 



z'=R(z+h) 



ar 



dz 



_ 1/ , l-2h 2h , l-4h 4A 
-R\ 1 + s 2 z 4 



d -l — # 

dz 



Then making use of the well-known relations 
*j) z m dz=0; if m^ — l 



and 



i£^=27ri;if m=-l 



it turns out, neglecting as usual terms containing powers 
of n, h, and a higher than the first, that 

P= P { -iP,== -ipfaT^l + (l W*>r 0 f(l + 6ft) 



or 



This last expression shows that, when the angle of attack 
is assumed small enough so that only the first power of a 
is retained, the component P{ of the lift vanishes in 
comparison with the component P,. 
Thus 



Lift=P„=AAr 



(34) 



This expression agrees in form with the corresponding 
one in incompressible flow with the auxiliary definition 

e. g., equation (27). 



CALCULATION OP THE MOMENT 

The moment arm QQ=m sin (<r— <p) and the force per 
unit length along the airfoil is pds (fig. 2). Hence the 
total moment about the origin 0 is given by 

M=<j)pmsin (<r— p)<fe=l(j)p(m sin a- cos <p 
— m cos o- sin <p)da 



Hence 



Now 



d£=ds cos a and drj= —ds sin a 
M= (^p(m cos <pd£+m sin^xftj) 
= <jW|+*fr?) =\ ^jprfm* 
^+^=i2.P.fc?f 



and since d£ =^r 



(35) 



Substituting for the pressure p the expression 
Co Mt .-I ft ( 1+ i,)(|*)' + I^(M) i + . . 

it follows that 

M=-i*(l+i„>.P.£(f)w 

Referring to the plane of the circle of radius It 

Performing the integrations in a manner analogous to 
that for the lift, it turns out that 

M=4xpoPW«[l-A+(l +§/^]- W?WM2 + Gft) 



(3G) 



or 



(37) 



This expression represents the moment about the origin 
of coordinates, and the moment about the center of the 
circle of radius R (into which the profile is mapped) can 
be obtained at once as 



1 | 7/j, 



or 



(38) 



where Mo=4Tp 0 *V i a , a is the moment about the center 
of the circle of radius R for the corresponding incom- 
pressible flow. 

If now d represents the distance of the center of 
pressure from the origin of coordinates, then 

M=Ld 

or 



d= 



4taW/i-/*+4^V) 

\ i L 



&a(l+3ktx) (39) 
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This expression shows that the airfoil has a constant 
center of pressure at a distance equal to 1/4 (1— Zh/i) of 
the chord from the leading edge. For a thin airfoil, 
say e=0.05 and for a stream velocity v 0 =0.8Z5 c a , the 
center of pressure, as compared with the corresponding 
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incompressible case, is nearer to the leading edge by 
about 2.5 percent of the chord. 
Langley Memorial Aeronautical Labobatoby, 
National Advisory Committee foe Aeronautics, 
Langley Field, Va., November 19, 1937. 



APPENDIX A 



APPLICATION OF THE THEORETICAL RESULTS 

As an example of the application of the theory to any 
particular case, the flows past a thin and a fairly thick 
symmetrical Joukowski profile for zero angle of attack 
will be calculated. Since no experimental results are 
available for purposes of comparison, it was considered 
hardly worth while to perform the rather lengthy and 
tedious calculations associated with angle of attack or 
circulation. 
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Fiotjbi 3.— Pressure and velocity distribution for the symmetrical Joukowski airfoil 
seottoo. «-=0.0S. n~0.7Q, o-O". 

Equation (23) for Avjv a simplifies considerably when 
the angle of attack a is taken to be zero. Thus 



A(12-2l A-4A 2 +30fe i -24A t -8A 6 ) . . 
gp ■ sm 5 



where 



_ (1+2A) (1-4^(1-/0' sing 



l—2k cos 8+k 1 
-h*k sin 5-A(cos S+h cos 2S) tan" 1 r^-—, - t 

1 — KC03 5 

-| (sin S+h sin 23) log (1-2A cos S+k l ) 



a* 

a. 
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Fiqcbe 4.— Pressure and velocity distributions for the symmetrical JoukowsM airfoil 
section e=0.ie. m"CU7, a=0*. 

The profiles chosen are defined by e=0.05 and «=0.15 
or k=l/21, k= 19/21; and A=3/23, k= 17/23, respect- 
ively. Tables I and II present the calculations and 
figures 3 and 4 show the velocity and pressure distri- 
butions for both incompressible and compressible flow. 
The values of /* chosen were the critical values ob- 
tained by plotting (v P lv 0 ) mtx against Void and then noting 
the intersection of this graph with that of (viv 9 )„ u 
against Pb/co as given by equation (9). Table III 
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presents the data and figure 5 shows the corresponding 
graphs. 

The expressions for Avfv a are given by: 



S 4 

























x 




















x 








































1 






















\ 








































































































I 

\ 
















-.05 


-.-./. 


5 




.< 


>35 




35 






Ji 














— 



.2 



.8 



1.0 



A £ 

Figure J. — Limiting values of ri/c, at zero angle ot attack for »-QJ0S and *~0.IS 

for =0.05 

^=^[0.05890 sin 5+0.05276 sin 25+0.00125 sin 35 
- 0.02585 J, (0) + 0.00246J S ] 



and 

for e=0.15 

— =ft[0^4905 sin 5+0.18053 sin 25+0.01151 sin 35 
— 0.40568,^(0) +0.10583lj] 

It is to be noted in tables I, II, and III that the maxi- 
mum velocity v P [vo for the incompressible flow occurs 
at about 5=35° and 5=45° for e=0.05 and e=0.15, 




o .1 .2 



.3 A .5 .6 .7 
(/chord 

€=.05 



8 .9 1.0 




t .2 .3 A 5 .6 .7 .8 J9 IJO 
(/chord 

Figukx 5 — Symmetrical Joukowskt profiles «~O.0J and c =-0.15. 

respectively. It is then assumed that the position of 
maximum velocity is independent of \l and maximum 
values for »p/p 0 are calculated for various values of ju. 
These values of (v F [vo)„ ax are given in table III and are 
used in obtaining the critical values of ft as shown in 
figure 5. The coordinates of the airfoils e=0.05, e=0.15 
are given in table IV and the corresponding contours 
in figure 6. 



APPENDIX B 



NOTATION 

t, y, rectangular coordinates in the plane of the 
circle. 

£, ri, rectangular coordinates in the plane of the 

profile. 
z=x+iy, r=$+iif 

r, 8, plane polar coordinates in the z plane. 

v x , v v , fluid velocity components along the x and y 
axes, respectively. 

i! e , tangential velocity on the circle. 

v P , tangential velocity on the profile correspond- 
ing to %. 

v= V^*+ V» magnitude of the fluid velocity. 
c, local velocity of sound in the fluid, 

p, density of the fluid. 

p, static pressure in the fluid. 

»oj <Jo> Po> Po, corresponding magnitudes in the undis- 
turbed stream. 

A», correction term to the velocity in incompres- 
sible flow due to compressibility. 
<t>, velocity potential of the incompressible flow. 

v r = — j£> component of velocity along the radius vector. 

Vt= — i ^> component of velocity perpendicular to the 
radius vector in the sense of 9 increasing. 

r 

R, radius of circle into which the profile is 

mapped. 
7, ratio of specific heats. 

w, complex potential of the incompressible flow 

in the f plane. 
a, angle of attack. 

210 



thickness coefficient of Joukowski profile (see 
%. 2). 



h= 



l+"« 



fc=i-p^=l-2A 

1 + 6 

r, "circulation about profile in the compressible 
fluid. 

r 0 , circulation about profile in the incompres- 
sible fluid. 

Ar=r— r 0 , contribution to the circulation due to com- 
pressibility. 
7, force vector on tlie airfoil. 
Pi, P„ components of P along the £ and 17 axes, 
respectively. 

M, moment about origin of coordinates in the 

plane of compressible flow. 
M e , moment about center of circle of radius i? 

in the plane of compressible flow. 
M 0 , moment corresponding to Af e in the plane of 

incompressible flow. 
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TABLE I 
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Case 1: e=(L05, ft=2y> fc = 2j"' «=0 
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Case 2: e =0.15, h=~, fc=g|. «=0 
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L350 


L870 


.6 


L143 






1.372 


I. 571 


.7 


1.151 






1.400 


1.366 


.8 


LI6I 






1. 430 


1.211 


LO 


L1S5 






L503 


L000 



l (±Y _41 + 2=l, wU h 7 .L408 
Kt/M ir+l/> T+l 



+0.16944 



TABLE IT 
COORDINATES OF THE AIRFOILS 





.-0.05 


•-0.15 


CdeV) 


i 


f 


i 






chord 


chord 


chord 


chord 


0 


i LOOM 


0.0000 


1.0000 


0.0000 


10 


.9930 


.0050 


.9938 


.0213 


20 


.9723 


.0155 


.9751 


.0411 


30 


.9471 


.0220 


.9444 


.0592 


40 


.8917 


.0270 


.9023 


.0739 


50 


.8339 


.0304 


.8195 


.0845 


60 


.7662 


.0318 


.7871 


.0907 


70 


.6901 


.0316 


.7163 


.0922 


80 


.6081 


.0295 


.6388 


.0891 


90 


.6228 


.0261 


.5558 


.0620 


100 


.4355 


.0218 


.4695 


.0715 


110 


.3496 


.0171 


.8830 


.0589 


120 


.2678 


.0124 


.2980 


.0151 


130 


.1923 


.0082 


.2178 


.0324 


140 


.1271 


.0049 


.1157 


.0211 


150 


.0732 


.0025 


.0840 


.0123 


160 


.0331 


.0010 


-03S8 


.0062 


170 


-006-1 


.0008 


.0093 


.0025 


180 


■-0000 


.0000 


.oooo 


.oooo 



i Kose. 



• Tall. 



